A Mutual Attraction Model for Both Assortative and Disassortative Weighted 

Networks 
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In most networks, the connection between a pair of nodes is the result of their mutual affinity 
and attachment. In this letter, we will propose a Mutual Attraction Model to characterize weighted 
evolving networks. By introducing the initial attractiveness A and the general mechanism of mutual 
attraction (controlled by parameter m), the model can naturally reproduce scale-free distributions 
of degree, weight and strength, as found in many real systems. Simulation results are in consistent 
with theoretical predictions. Interestingly, we also obtain nontrivial clustering coefficient C and 
tunable degree assortativity r, depending on m and A. Our weighted model appears as the first one 
that unifies the characterization of both assortative and disassortative weighted networks. 
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The past few years have witnessed a great deal of in- 
terest from physics community to understand and char- 
acterize the underlying mechanisms that govern complex 
networks. Prototypical examples cover as diverse as the 
Internet Q, the World-Wide Web 0], the scientific col- 
laboration networks (SCN) and world-wide airport 
networks (WAN) 0- Q As a landmark, Barabasi and Al- 
bert (BA) proposed their seminal model that introduces 
the linear preferential linking to mimic the topological 
evolution of complex networks [7]]. However, networks 
are far from boolean structure. The purely topological 
characterization will miss important attributes often en- 
countered in real systems. For example, the amount of 
traffic characterizing the connections of communication 
systems or large transport infrastructure is fundamental 
for a full description of these networks @ . This thus calls 
for the use of weighted network representation, which is 
often denoted by a weighted adjacency matrix with el- 
ement Wij represents the weight on the edge connect- 
ing vertices i and j. In the case of undirected graphs, 
weights are symmetric Wij = Wji , as this letter will focus 
on. A natural generalization of connectivity in the case 
of weighted networks is the vertex strength described as 
Si = J2j£T(i) w ij> wnere the sum runs over the set T{i) 
of neighbors of node i. This quantity is a natural mea- 
sure of the importance or centrality of a vertex in the 
network. Most recently, the access to more complete em- 
pirical data and higher computation capability has al- 
lowed scientists to consider the variation of the connec- 
tion weights of many real graphs. As confirmed by mea- 
surements, complex networks not only exhibit a scale- free 
degree distribution P(k) ~ k~i with 2< 7 <3 H@, but 
also the power-law weight distribution P(w) ~ w~° 
and the strength distribution P(s) ~ s~ a Highly 
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correlated with the degree, the strength usually displays 
scale- free property s ~ k$ with (3 > 1 ja, [nj, uM ■ Moti- 
vated by those findings, Alain Barrat et al. presented a 
model (BBV for short) to study the growth of weighted 
networks [13 ■ Controlled by a single parameter 5, BBV 
model can produce scale-free properties of degree, weight 
and strength. But its disassortative property [l^, uM 
(i.e. the hubs are primarily connected to less connected 
nodes), as observed in real technological and biological 
networks, can hardly give satisfying interpretations to 
social networks like the SCN, where the hubs are very 
likely to be linked together (assortative mixing). Previ- 
ous models as far as our knowledge can generate either 
assortative networks 0, 0, 0] or disassortative ones 
|l2l fl3l FLU ITU], but rarely both. Thus, some questions 
arises here: why are social networks all assortative, while 
all biological and technological networks opposite? Is 
there a generic explanation for the observed incompati- 
ble patterns, or does it represents a feature that needs 
to be addressed in each network individually? Our work 
may shed some new light to these questions. 

Former network models often impress on us such a net- 
work evolution picture: pre-existing nodes are passively 
attached by newly added node according to the preferen- 
tial linking mechanism. This scenario, however, lacks the 
other side of fact that old nodes will choose the young 
at the same time. In addition, this evolution picture also 
ignores the universal mutual attraction between existing 
components, which leads to the creation and reinforce- 
ment of connections. This idea has been partly reflected 
in the studies of Dorogovtsev and Mendes (DM) who 
proposed a class of undirected and unweighted models 
where new edges are added between old sites and exist- 
ing edges can be removed. In this letter, we will present 
a model for weighted evolving networks that considers 
the topological evolution under the general mechanism of 
mutual attraction between nodes. In contrast with previ- 
ous models where weights are assigned statically 
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or rearranged locally [T2I Il3| , our model allow weights 
to be widely updated. It can mimic the reinforcement 
and creation of internal links as well as the evolution 
of many infrastructure networks. Specifically, the model 
can generate a diversity of scale- free quantities, nontrivial 
clustering property, and tunable assortativity coefficient. 
Therefore, one can easily find explanations to various real 
networks by our microscopic mechanisms. 

The model starts from an initial configuration of No = 
m isolated nodes with no connections between each other. 
At each time step, a new isolated node n is introduced 
into the system. Then every existing node i (including 
the newly-added one) selects m other existing nodes for 
potential interaction with the probability 

« -• A 
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To guarantee the isolated new node be chosen by oth- 
ers the model requires A > 0, and hence A is called the 
initial attractiveness |2l| . If two unconnected nodes are 
mutually selected, then an internal connection is created 
between them. If there already exists a connection be- 
tween them, their link is just strengthened by increasing 
weight wq — 1. Here, m is the number of candidate nodes 
(of each site per step) to create and strengthen connec- 
tions, and the initial attractiveness A governs the prob- 
ability for "young" nodes to get new links and weights. 
After the weights have been updated, the growth process 
is iterated by introducing a new vertex until the desired 
size of the network is reached. 

One can easily find appropriate interpretations to real 
networks from our model mechanism. Take the SCN 
for example: the collaboration of scientists require their 
common interest and mutual acknowledgements. Unilat- 
eral effort does not promise effective activity. Speaking 
by the model language, though the nodes with low degree 
would like to connect to nodes with large degree, the lat- 
ter do not necessarily wish to be linked by former. On the 
other hand, two scientists both with strong scientific po- 
tentials (large strengths) and long collaborating history 
are more likely to publish papers together. The above 
description of our model could also satisfactorily explain 
the WAN where the edge weigh denotes the relative mag- 
nitude of the traffic along a flight line. During the evolu- 
tion of WAN, the airlines are more likely to open between 
metropolises that hold a high status in both economy and 
politics (with large strengths). With the improvement of 
economy and the expansion of population, the air traffic 
between connected metropolises will increase much faster 
than that between smaller cities. Due to their impor- 
tance, there is an obvious need for smaller cities to build 
new airports to connect those metropolises. But the limit 
of energy and resources leads to the fact each node can 
only afford a finite number of connections. Therefore, 
they have to choose in front of the vertex pool. 

The model time is measured with respect to the num- 
ber of nodes added to the graph, i.e. t — N — No, and 
the natural time scale of the model dynamics is the net- 
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FIG. 1: Numerical simulations by choosing A = 1. Data are 
averaged over 10 independent runs of network size N = 8000: 
(a) Cumulative probability strength distribution P(s) with 
various values of m. Data are consistent with a power-law 
behavior P(s) ~ s~ a . The inset reports the values of a ob- 
tained by data fitting (full circles) in comparison with the 
theoretical prediction a = 2 + A/m 2 (line), (b) Cumulative 
probability degree distribution P(k) with m — 1 and m = 2. 
Data fitting confirms its scale-free property, (c) Cumulative 
probability distribution of weight with different m, in agree- 
ment with the power-law tail P(w) ~ w~ e . As shown in its 
inset, the data fitting also gives values of (full circles) as 
predicted by analytical calculation (line), (d) The average 
strength Sj of nodes with connectivity kt for different m. In 
the log-log scale, we observe the nontrivial strength-degree 
correlation s ~ k 13 , with the exponent j3 versus network size 
N (see the inset). 



work size N. Using the continuous approximation, we 
can treat k,w,s and the time t as continuous variables 
0, 0| ■ The time evolution of the weights toy can be 
computed analytically as follows: 



dwij 
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m 2 ( Sl + A)(sj + A) 



Hence, the strength Sj(i) is updated by this rate: 
dsi ^ dwij m 2 (si + A) m 2 (si + A) 



dt 



(2) 



(3) 
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The last expression is recovered by noticing that 

J2(si+A) = V Si+At = [ ^i^dt+At = (m 2 +A)t. 
i Jo at 

From Eq. (3), one can obtain the scaling of Si(t) versus 
t as Sj(t) ~ t x , which also implies the scale-free distri- 
bution of strength P(s) ~ s~ a with the exponent [l^ 
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FIG. 2: Numerical simulations by choosing A = 5. Data 
are averaged over 10 independent runs of network size JV = 
8000: (a) Cumulative probability strength distribution P(s) 
with various values of m. Data are consistent with a power- 
law behavior P(s) ~ s -c \ The inset reports the values of 
a obtained by data fitting (full circles) in comparison with 
the theoretical prediction(line) . (b) Cumulative probability 
distribution of degree for different m. Data fitting confirms its 
scale-free property: P(k) ~ k (c) Cumulative probability 
distribution of weight with different m, in agreement with 
the power-law tail P(w) ~ w~ e . As shown in its inset, the 
data fitting also gives values of 8 (full circles) as predicted 
analytically (line) . (d) The average strength s, of nodes with 
connectivity ki for different m. In the log-log scale, we observe 
the nontrivial strength-degree correlation s ~ k 13 , with the 
exponent (5 versus network size N (see the inset). 
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One can also obtain the evolution behaviors of weight 
and degree, and hence their power-law distributions |22| : 
P(w) ~ w~ e with 



2A 



(5) 
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and P(fc) ~ fc -7 with 7^2 + A/to 2 = a as < — > oo. 

We performed numerical simulations of networks gen- 
erated by choosing different values of A and to and recov- 
ered the theoretical predictions. We have also checked 
that the scale-free properties of our model networks 
are almost independent of the initial conditions. Fig. 
l(a)-(d) report the probability distributions of strength, 
weight and degree, as well as the strength-degree correla- 
tion, fixed A = 1 and tuned by to. Specifically, Fig. 1(a) 
gives the probability distribution P(s) ~ s~ Q , which is in 
good agreement with the theoretical expression [Eq. 4]. 
Probability weight distribution also recovers the power- 
law behavior P(w) ~ w~ e [Fig. 1(b)] with 9 as predicted 



FIG. 3: (a) Clustering coefficient C depending on both m 
and A with network size iV = 8000. (b) Degree assortativity 
r depending on both m and A with network size iV = 8000. 



analytically [Eq. 5] . Fig. 1 (c) shows the scale- free degree 
distribution P(k) ~ fc~ 7 and Fig. 1(d) reports the av- 
erage strength of vertices with degree ki, which displays 
a nontrivial power-law behavior s ~ fc' 3 as confirmed by 
empirical measurements. The inset of Fig. 1(d) indicates 
that the exponent f3 decreases slowly with the network 
size [2^ , which is noticeably different from the linear cor- 
relation (P = 1) as obtained in most previous models. 
Again, Fig. 2(a)-(d) show the simulation results by fixing 
A = 5 and adjusting m. In comparison with Fig. l(a)(b), 
the distributions of strength and degree for A — 5 both 
behave exponential corrections in the zone of low degree. 
This phenomenon occurs at large A and the exponential 
part are very similar with the empirical findings in some 
social networks like SCN |2^|. In the zone of large de- 
gree, however, we can still observe the scale-free behavior 
which again recover the theoretical exponent expressions. 
It is worth remarking ahead that the model at large A 
can generate the assortative property too, which is seen 
in social networks. Thus, the introduction of A is impor- 
tant for our model to mimic social networks. 

To better understand the degree correlations of our 
model networks, we also studied the clustering coeffi- 



cient C (which describes the statistic density of connected 
triples) |2i| and degree assortativity r |24| depending on 
the model parameters A and m. As presented in Fig. 
3(a), C for fixed m monotonously decreases with A, and 
C for fixed A monotonously increases with m. Generally, 
it can be adjusted in the range ]0,1[. Obviously, the clus- 
tering property of our model is tunable in a broad range 
by simultaneously varying m and A, which makes it more 
powerful in modelling real networks. As shown in Fig. 
3(b), degree assortativity r for fixed m, unlike the clus- 
tering case, increases with increasing A; while r for given 
A decreases with m. For small A and large m, the model 
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generates disassortative networks which can best mimic 
technological networks like the Internet [ij and WAN or 
even biological networks. While at large A and small 
m, assortive weighted networks emerge and can be used 
to model social graphs as the SCN. Actually, enhancing 
the initial attractiveness A will considerably increase the 
chances for "young" nodes to be linked and strength- 
ened. As low-degree nodes take the majority in the sys- 
tem, larger A will lead to the stronger affinity between 
"young" vertices, and thus they can link together more 
easily. This explains the origin of assortative mixing in 
our model and may also shed some light on the old open 
question: why social networks are different from other 
networks in degree assortativity? Considering humans 
are active elements, it is plausible that the components of 
social networks possess considerable initial attractiveness 
(large A). On the other end, as m controls the interac- 
tion level of internal connections, increasing m will make 
the hubs become busier and busier, as they have to be 
linked by more and more "young" sites. It may explain 
why the disassortativity of the model is increasingly sen- 
sitive to m. Combining these two parameters together, 
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the mutual attraction model integrates two competitive 
ingredients that may be responsible for the mixing dif- 
ference in complex networks. 

The universal mutual attraction between nodes and 
the existence of initial node attractiveness are two im- 
portant united ingredients of our model in mimicking real 
weighted networks. The general dynamics of node inter- 
action proposed in this letter provides a wide variety of 
scale-free behaviors, nontrivial clustering coefficient and 
tunable degree assortativity. As far as our knowledge, 
the weighted network model presented here appears as 
the first one that can both mimic assortative and disas- 
sortative networks under a unified evolution dynamics. 
Its obvious simplicity and reproduced real-world variety 
allow more specific mechanisms to be integrated into fu- 
ture modelling efforts. Above all, the Mutual Attrac- 
tion Model we presented here implies us the possible and 
worthwhile efforts in exploring the unified mechanisms 
behind various networks. 
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